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ABSTRACT 



Context. Comet 17P/Holmes underwent a dramatic outburst in October 2007, caused by the sudden fragmentation of 
its nucleus and the production of a large quantity of grains scattering sunlight. 

Aims. We report on 90 GHz continuum observations carried out with the IRAM Plateau de Bure interferometer on 27.1 
and 28.2 October 2007 UT, i.e., 4-5 days after the outburst. These observations probed the thermal radiation of large 
dust particles, and therefore provide the best constraints on the mass in the ejecta debris. 

Methods. The thermal emission of the debris was modelled and coupled to a time-dependent description of their 
expansion after the outburst. The analysis was performed in the Fourier plane. Visibilities were computed for the two 
observing dates and compared to the data to measure their velocity and mass. Optical data and 250-GHz continuum 
measurements published in the literature were used to further constrain the dust kinematics and size distribution. 
Results. Two distinct dust components in terms of kinematic properties are identified in the data. The large-velocity 
component, with typical velocities Vo of 50-100 m s" 1 for 1 mm particles, displays a steep size distribution with a size 
index estimated to q = -3.7 (±0.1), assuming a minimum grain size of 0.1 /im. It corresponds to the fast expanding shell 
observed in optical images. The slowly-moving "core" component (Vo = 7-9 m s _1 ) detected near the nucleus has a size 
index \q\ < 3.4 and contains a higher proportion of large particles than the shell. The dust mass in the core is in the 
range 0.1-1 that of the shell. Using optical constants pertaining to porous grains (50% porosity) made of astronomical 
silicates mixed with water ice (48% in mass), the total dust mass M^ust injected by the outburst is estimated to 4-14 
x 10 11 kg, corresponding to 3-9% the nucleus mass. 

Key words. Comet: individual: 17P/Holmes - Radio continuum: solar system - Techniques: interferometric 



1. Introduction 

Comet 17P/Holmes is a periodic comet of the Jupiter fam- 
ily that orbits the Sun with a period P = 6.9 years. It 
passed perihelion on 4 May 2007 at 2.05 AU from the Sun. 
On 24 October 2007, at 2.44 AU from the Sun and 1.63 AU 
from the Earth, the comet suddenly increased in bright- 
ness from a total visual magnitude m v ~ 17 to 2.5 to be- 
come a naked-eye object for months (lGreenll2007l ; ISekaninal 
2009). Comet 17P/Holmes underwent a similar outburst 
shortly before 6 November 1892, at the time of its discov- 
ery, followed by a similar e vent on 16 January 1893 (see 
the review of ISekaninall2009D . After these events, the comet 
appearance at large scales was a bubble-lik e shape quickly 
expanding into interpl anetary space (e.g., iMontalto et al.l 
l2008tlHsieh et al.ll2010D . The onset of the 2007 outburst oc- 
curred probably near 23.3 October UT (Hsi eh et al.ll2010h 
with the peak of o ptical brighteni ng observed around 25.0 
October UT (e.g., iLi et all 



tering sunlight. Determining the amount of material that 
split off from the nucleus and the size distribution of the 
particle debris in the cloud of dust ejecta is important to 
constrain the origin of the fragmentation process. Though 
sizeable indivi dual fragments radiating outwards were pos - 
sibly observed ([Gaillard et al.|[2"007t IStevenson et al.l l2010). 
the huge cross-sectional area of dust scattering sunlight sug- 
gests that the dust in the comet Holmes 2007 outburst was 
dominated by small particles. 

The potential of millimetre and submillimetre- 
wavelength continuum observations for the study of 



20111) . These outbursts were 



likely caused by a sudden fragmentation of the nucleus, fol- 
lowed by the production of a large quantity of grains scat- 



* Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 



come t ary dust has been demonstrated (e.g.. iJewitt fc Luul 
1990 . 119921: 1 Jewitt k MatthewsHl999tlAltenhoff et al.l ll999: 
iBockelee-Morvan et al.l l2010bj ). Such observations probe 
the thermal radiation of millimetre-sized dust particles, 
and therefore usefully complement optical and infrared 
observations, that are sensitive to micrometric particles. 
By measuring the radiation from large particles, they are 
of high value to measure dust masses. Constraints on the 
dust properties, e.g., the size distribution, can be obtained 
if measureme nts of the spectral index o f the dust emission 
are available (jJewitt fc Luulll990l[l992l ). 

Single-dish continuum measurements at 1.1 mm wave- 
length of comet 17P/Holmes obtained with the 30-m tele- 
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scope of the Institut de Radi oastronomie mil li met riquc 
(IRAM) have been presented bv lAltenhoff et al.l <|2009h . A 
short report of observations at 1.3 and 0. 8 mm conducte d 
with the Submillimeter Array was given bv lQi et all ([20101) . 
In this paper, we present interfcromctric continuum ob- 
servations performed at 3.3 mm wavelength on 27 and 28 
October 2007 UT. These observations, carried out with the 
IRAM Plateau de Bure interferometer, provided images of 
the dust coma at ~ 6" angular resolution which corresponds 
to 7100 km diameter at the comet distance. The observa- 
tions and data products are presented in Sects. [2] and [3] 
They were analysed with a model of dust thermal emission 
coupled with a time-dependent model of the expansion of 
the cloud. The modelling approach is described in Sect. 0] 
and the data analysis is done in SectJSJ A discussion of 
the results obtained on the properties of the dust ejecta 
follows in ScctlHl A prelimi n ary r e port of the observations 
was given bv lBoissier eTail (|2008l 200j|). 

2. Observations 

Observations were undertaken on 26.92-27.29 and 28.08- 
28.33 October 2007 UT at the IRAM Plateau de 
Bure interferometer (P dBI) situated in the French Alps 
([Guilloteau et al.l [19921. Therefore, they were performed 
soon afte r the annou ncement of the comet outburst on 24 
October (lGreenll2007ft as a Target of Opportunity program. 
The comet was at a geocentric distance A ~ 1.63 AU and 
a heliocentric distance r/j ~ 2.45 AU. The Sun was at a po- 
sition angle p. a. = 38-39° on the plane of the sky, and the 
phase (Sun-Comet-Earth) angle was ~ 16°. The comet was 
tracked using the ephemeris provided by the HORIZONS 
system (solution JPL K077/6). At the time of the obser- 
vations, the six 15-m antennas of the interferometer were 
set in the compact D configuration. On 27 October, only 
five antennas were available. The baseline lengths projected 
onto the plane of the sky ranged from 15 to 100 m. The 
UD-coverage acquired over the course of observations is pre- 
sented in Fig. [T] 

The observations were performed using the 3-mm dual 
polarization (4-GHz bandwidth) receivers tuned to the fre- 
quencies of the J(l-0) HCN (88.6 GHz) and J(l-0) HNC 
(90.6 GHz) lines on 27 and 28 October, respectively. The 
two polarizations were observed using the 8 units of the 
Plateau de Bure correlator: for each polarization, one unit 
of 20 MHz bandwidth was centred on the molecular line, 
and three 320 MHz units were placed nearby for measur- 
ing the continuum emission. The effective total bandwidth 
for the continuum observations was thus about 2 x 0.9 
GHz. The observing cycle was typically: pointing, focussing, 
cross-correlation on a calibration source, 2 min of auto- 
correlation (ON-OFF measurements, 1 min on source, for 
line observations with the angular resolution of the primary 
beam) and 30 scans (45 s each) of cross-correlation on the 
comet. A log of the observations is presented in Tabled] The 
angular resolution of the synthesized beam is ~ 6.4" and 
5.5" on 27 and 28 October, respectively and the Half Power 
Beam Width (HPBW) of the primary beam is about 54" 
at the observed frequencies. The analysis of the HCN and 
HNC data (cross-correlation, i.e. interferometric mode, and 
ON-OFF measurements) will be presented in a forthcom- 
ing paper. 

The data r eduction w as done using the GILDAS soft- 
ware package (|Petvi I2005D . For both datasets (27 and 28 
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Fig. 1. uw-coverage for 27 and 28 October. The (u,v) 
points in black represent the part common to the two dates. 



October) the bandpass calibration was carried out observ- 
ing 3C454.3. The instrumental phase and amplitude vari- 
ations were deduced from observations of the calibrator 
0355+508. The weather conditions at the Plateau de Bure 
were outstanding, and the resulting phase rms is at most 
20°, with values lower than 15° for most of the antenna 
pairs. Finally the absolute flux density scale was determined 
from measurements of 3C84 and MWC349 fluxes. The con- 
sistency (within 2.5%) of the calibrator fluxes measured on 

27 and 28 October implies an excellent relative calibration 
between the two dates. The uncertainty in the absolute flux 
calibration is about 15% for both days. 

The interferometric maps deduced from the whole data 
set are presented in Fig. [5] and their main characteristics 
(flux density at the centre of the map and astrometric po- 
sition of the peak of brightness) are reported in Tables [T] 
and [5] Note that the maps have been centred on the posi- 
tion of the peak of brightness. The comparison of the 27 and 

28 October maps is difficult due to differences in the shape 
and size of the interferometric beam, caused by different uv- 
coverages. Maps obtained with a similar ww-coverage using 
restricted data sets (from 2h00 to 7h00 UT, same 5 anten- 
nas for both dates) are also presented in Fig. [21 and their 
characteristics are given in the second part of Tables [T]-[2] 



3. Data analysis 

The detected 3.3 mm emission is due to the thermal emis- 
sion of dust particles in the coma. Indeed, a flux density 
of 2 mJy at 90 GHz would correspond to a photometric 
diameter of 60 km for a spherical slow rotator at the equi- 
librium temperature of 207 K expected at 2.45 AU from 
the Sun. Instead, the effective diameter of the nucleus of 
17P / Holmes is estimated to be 3.2 km ([Snodgrass et al.l 
I2006h . 

The first interesting feature in these 3.3-mm continuum 
measurements is the small variation of the flux density in 
24 hours elapsed time. The peak flux density within the ~ 
5-6" field-of-view decreased by 20-25% only between the 
two dates (Table [J). 



3.1. Astrometry and azimuthal variations 

As seen in Table [5J the astrometric position of the max- 
imum brightness (O) almost coincides with the position 
of the nucleus given by the ephemeris (C). An (O-C) off- 
set of ~ +0.3" at the l-cr confidence level is measured 
in Declination for both dates when considering the whole 
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Table 1. Observation log and flux measurements. 





Date UT 
(October 2007) 


Freq. 
(GHz) 


N A a 


Beam 
(" x ") 


S/T A b 

(Jy/K) 


Peak intensity 
(mjy/beam) 


S/N 


AO c 
(") 


Total datasets 

26.92-27.29 

28.08-28.33 


88.6 
90.6 


5 
6 


6.88x5.97 
6.44x4.77 


23.14 
22.4 


2.38 ± 0.07 
1.9 ±0.1 


34 
19 


0.2 
0.3 


Subsets e 

27.08-27.29 

28.08-28.29 


88.6 
90.6 


5 
5 


7.49x5.89 
7.19x5.66 


23.14 
22.4 


2.45 ± 0.12 
2.06 ±0.14 


20 
15 


0.3 
0.4 



a Number of observing antennas. 

b Conversion factor between flux and antenna temperature. 

c AO is the astrometric precision of the position measured on the map. It is given by the interferometric beam divided by the 

signal to noise ratio. 
d In this part we take into account all the observations for each day. 

e Here the maps were done with the observations performed in the common UT range of 27 and 28 October. 



Table 2. Position of the brightness centre. 





Date 


Freq. 


UT time" 


Ephemeris Position (C) 6 ' c 


Observed Position (O) 


(0-C) d 








RA 


Dec 


RA 


Dec 


(RA, Dec) 


(2007) 


(GHz) 


00 


(h:min:s) 


(«:':") 


(h:min:s) 


(":':") 


(",") 


Total datasets 
















27 October 


88.6 


7.0000 


03:51:30.418 


50:19:16.48 


03:51:30.418 


50:19:16.76 


(±0.00,±0.28) 


28 October 


90.6 


8.0000 


03:50:30.701 


50:23:04.37 


03:50:30.733 


50:23:04.66 


(±0.31,±0.29) 


Subsets 
















27 October 


88.6 


7.0000 


03:51:30.418 


50:19:16.48 


03:51:30.436 


50:19:16.59 


(±0.17,±0.11) 


28 October 


90.6 


6.0000 


03:50:35.552 


50:22:46.94 


03:50:35.569 


50:22:47.26 


(±0.16,±0.31) 



a Reference time for the position in RA, Dec of the comet. 

6 The coordinate system is apparent positions, geocentric frame. 

c Ephemeris computed from the HORIZONS system, solution JPL K077/15 based on 3265 astrometric data from 1964 July 16 
to 2008 April 22. 

6 Offset between the peak of the continuum emission and the position given by the ephemeris. 



data set. On optical images, the inner coma of 17P/Holmes 
presents an asymmetric distribution, mainly characterized 
by the presence at the position angle p. a. = 220° (cor- 
responding to the anti-Sun direction) of a bright elon- 
gated cloud (so-called "blob"), quickly separating at a 
rate of 9-10 "/day from the (much brig hter) condensa- 
tion of material surroundin g the nucleus flMontalto et al.l 
[2^f31 IWatanabe et al J 120091; iReach et al.ll2010D . The (O-C) 
marginal offsets measured in the PdBI maps (p. a. within 0- 
60°) are not towards the direction of this expanding blob: 
besides the uncertainty of the measurements and of the 
comet ephemeris, they can be due to spatial asymmetries 
of the brightness distribution in either direction. At the 
time of our observations, the blob was at a projected dis- 
tance of - 1 7" (27 October) and 27" (28 October) from 
the nucleus (jMontalto et alj l2008t IWatanabe et al.ll2009j) . 
i.e., at the edge of the primary beam of the PdBI anten- 
nas that defines the extent of the interferometric map. Its 
thermal emission is not seen in our images (Fig. [5]). Thus, 
the amount of millimetre-sized dust particles in this blob 
was too low with respect to the amount of material sur- 
rounding the nucleus for a detection in the PdBI maps. 
We have also to consider that more massive particles have 
lower expansion velocities and are less subject to radiation 



pressure, so that this bright blob of micrometric particles 
seen in optical and infrared images was likely deficient in 
particles radia ting at 3.3 mm. C omparing optical and in- 
frared images, IReach et al.1 (|2010f ) concluded that the blob 
comprised particles of intermediate (~ 10 /zm) sizes. 

Spatial asymmetries are marginally present in the PdBI 
maps. To enhance spatial features, we subtracted a sym- 
metric image to the measurements. The subtraction was 
done in the Fourier plane, that is, we subtracted the vis- 
ibilities characterizing the symmetric image to the mea- 
sured visibilities. Interferometric maps were then produced 
with the new sets of visibilities (i.e., uw-tables). The vis- 
ibility amplitude V of the symmetric images follows V = 
0.052r-° 83 Jy and V = 0.025r.-°- 66 Jy, for 27 and 28 
October, respectively. This corresponds to the angular av- 
erage of the data in the Fourier plane described in Sect. 13.21 
Figures|31i, h show, for the two observation dates, the resid- 
ual interferometric images. We also plot residual images 
obtained by applying a factor / = 0.8 (Fig. |3p, f) and 
0.9 (Figs. [5g) to visibilities of the symmetric images. 
An emission feature is observed at the 2a level South-East 
from the nucleus position (p. a. ~ 160°), in the 27 October 
residual images obtained for / = 0.9 and 1.0. As a matter of 
fact, this feature is discernable on the 27 October original 
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Fig. 2. Intcrfcromctric maps of the 3.3 mm continuum 
emission of comet 17P/Holmes obtained on 27 and 28 
October 2007 UT with PdBI. Top: full set of measurements. 
Bottom: maps constructed using the same uv-coverage, i.e., 
from data acquired with 5 antennas and between 2h00 and 
7h00 UT for both dates. The interferometric beam is plot- 
ted in the bottom left corner. The maps have been centred 
on the position of maximum brightness. The arrow rep- 
resents the Sun direction projected on the sky. The level 
spacing is 2a, with a — 0.093 (top figure) and 0.10 (bot- 
tom) mJy/beam for 27 October, and a = 0.12 (top) and 
0.14 (bottom) mJy/beam for 28 October. 



map (Figs. [2] left and[3^i). Its intensity is ~ 10% the peak 
intensity in the centre of the 27 October map. This feature 
is not present in the 28 October images. More marginal 
features at the noise level (la) are observed at position an- 
gles corresponding approximatively to the Sun direction (27 
October), and to the tail direction (28 October) (p. a. (Sun) 
= 38°, p. a. (tail) = 218°). All these features are also ob- 
served when the restricted data sets are used. We show in 
Sect. [5] that the velocity of the 1-mm sized particles in the 
expanding shell is 50-100 m s -1 . This is compatible with 
the disappearance of the southward feature in the PdBI 
map of 28 October. 

3.2. Radial distribution 

The data provide information on the radial distribution of 
the emission at projected distances D of typically 3500 to 
20 000 km from the nucleus (this corresponds to the range 
of X/Bi projected at the distance of the comet, where Bi 
is the baseline length and A is the wavelength). For comets 
observed in steady-state activity, the column density and 
brightness distribution of the dust in the inner coma is ex- 
pected to vary in first approximation according to 1/D, 
not considering deviations from this ideal law caused by ra- 
diation pressure, particle fragmentation, and asymmetries 
related to the nucleus outgassing geometry. There are sev- 
eral examples of surface brightness profiles at millimetre 
and submillimetre wavelengths consistent with this ideal 
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Fig. 3. Asymmetries in the interferometric 3.3 mm maps 
of comet 17P/Holmes from image processing. Results for 
27 and 28 October are plotted in left (a— d) and right (e— 
h) panels, respectively. Top figures (a, e): observed images 
as in Fig. [5J All other figures have had a symmetric im- 
age subtracted (in the Fourier plane) to enhance azimuthal 
asymmetries. The symmetric images for 27 and 28 October 
are described by the fits of the visibilities given in Fig. @] 
multiplied by a factor / equal to 0.8 (b, f), 0.9 (c, g), and 
1.0 (d, h). The interferometric beam is plotted in the bot- 
tom left corner. The level spacing is lcr, except for the top 
figures (2a, see Fig. [2]), with a — 0.093 mJy/beam for im- 
ages (a— d), and a — 0.12 mJy/beam for (e— h). Dashed 
contours correspond to negative fluxes. 



law (e .g. JJewitt fc Matthewslll999tlBockelee-Morvan et al.l 
l2010bl ). On the other hand, when observations are con- 
ducted soon after a massive outburst, the radial profile may 
strongly differ from this law. The profile may resemble that 
of a point-like source, if the size of the expanding cloud is 
much smaller than the angular resolution of the interfero- 



4 



J. Boissier et al.: The 3.3-mm continuum emission of comet 17P/Holmes after its 2007 outburst 

Table 3. Optical constants m = n — ik at 3.3 mm. 



10 



10 




Compound 


n 


k 




Reference 


Water ice 


1.79 


2.00x10" 


-3 


Warren & Brandt (2008) 


Astro, silicates 


3.44 


1.50x10" 


-2 


Draine (1985) 


Forsterite 


2.63 


2.06x10" 


-4 


Fabian et al. (2001) 


Organics 


2.28 


2.50x10" 


-3 


Pollack et al. (1994) 


Silicate mixt. a 


2.76 


2.82x10" 


-3 


this work 



20% astronomical silicates and 80% forsterite. 



-radius [m] 



visibility at r uv 
ties. 



> 60 m is constant, within the uncertain- 



Fig. 4. Visibilities measured on 27 (black symbols) and 
28 (red symbols) October. Dashed lines show fitted de- 
pendences with w-radius: 0.052r~° 83 (27 October) and 
0.025r ra - 66 (28 October) in janskys. 



metric beam. Taking into account the elapsed time between 
the outburst onset time and the PdBI observations, this 
would happen for grains with velocities significantly less 
than 10 m/s. Conversely, if the millimetre-sized particles 
that contribute to the emission have sufficiently high ve- 
locities, a ring-like structure, arcs or bright condensations 
(not coinciding with the nucleus position) can be observed, 
depending on the geometry of the expansion. In the present 
case, this is not observed. Finally, a variety of profiles, even- 
tually approaching the l/D law in restricted regions of the 
coma, can be observed if the outburst was followed by a sus- 
tained production of millimetre-sized particles (e.g., due to 
dust fragmentation) and/or if the particles contributing to 
the signal have a broad range of velocities. The decline of 
the dust emission from 27 and 28 October is governed both 
by the production curve of the particles and their velocity. 

As discussed in iBoissier et al.l (|2007|) . an appropriate 
way to study the radial distribution of material observed 
by intcrfcromctric techniques is from the analysis of the 
visibilities which are the direct output of interferometers. 
This method allows us to avoid uncertainties in the maps 
resulting from the partial wii-coverage (see Fig. [lj . In ad- 
dition, due to filtering and missing data at short antenna- 
spacings, maps lose information concerning the comet emis- 
sion at large scales (see Fig. [2] where the emission is barely 
detected at D > 10"). Observed visibilities have been az- 
imuthally averaged in the uw-plane to investigate the ra- 
dial brightness distribution. Figure [4] displays the visibility 
amplitudes V as a function of w-radius r uv (the baseline 
length projected onto the plane on the sky) for the two 
dates. They vary according to r-°- 83±0 06 and r-^ 66±0 07 
for 27 and 28 October, respectively. The radial variation of 
the visibilities is indicative of a more compact distribution 
than expected for a l/D variation of the column density. 
For such a distribution, a r~} variation would have been 
observed (Bockelcc-Morvan et al. 2010a) (V independent of 
Mf-radius characterizes compact sources unresolved by the 
interferometic beam). A 23% decline of the emission is ob- 
served at r uv < 30 m between the two dates, whereas the 



4. Modelling 

4.1. Dust thermal emission 

Models of the thermal emission of cometar y dust in the 
microwave domain have been presented by Ijewitt fc Luul 
(|1990LI1992D . The radiation of dust grains at A = 3.3 mm de- 
pends of the grain absorption efficiency factor Q a b s , which 
is a function of the grain size and of the complex refrac- 
tive index m\ of the material. We calculated Q a bs at A 
= 3.3 mm for different mixtures of silicates, organics and 
water ice using the Mie theory for spherical and homoge- 
nous grains. Infrared signatures of water ice grains were 
observed in the coma of 17P/Holmcs, and lifetime calcula- 
tions indicate that dirty and pure ice grains with a >100 
/ im survived for several days or more after the outburst 
(jYang et al.l l2009h . We also considered the porosity of the 
grains. The effective m\ of fluffy grains was computed in 
a two-step process. First, we computed the refractive in- 
dex of basic units made of silicates and water ice using 
the M axwell Garnett effective m edium theory (EMT) fol- 
lowing [Greenberg^^HaSa (|1990f ). Mixtures of organics and 
ice, as well as mixtures of different silicates were also con- 
sidered (Table [3]). The effective refractive index of the fluffy 
particles was then deduced using the Maxwell Garnett for- 
mula for two-component mixtures, taking vacuum for the 
core material and the silicate/ice or organic/ice mixture in 
the matrix (hollow sphere). Thi s prescription was used by 
iKrueeel fc Siebenmorgenl (|l994h for modelling Q a b s of in- 
terstellar dust. EMT theories are only an approximation to 
the real optical behavior of composite media, and give more 
accurate results when the size of the grains is small relative 
to the wavelength, and when the volume of the inclusions 
(here silicat es, organics and voids) is small with r espect to 
the matrix (JOsscnkoof 199ll lPerrin fc Lamvlll996] ). The use 
of more exact theories is beyond the scope of this paper. 

The refractive indices used in this study are given in 
Tablc[3] We considered different materials, with optical con- 
stants which might be re presentative of cometary material 
(|Hanner fc Bradlevll2005D . Optical constants at 3.3 mm are 
rare in the literature. The value for the silicate mixture 
made of forsterite and astronomical silicates was computed 
with the Maxwell Garnett formula, with forsterite consti- 
tuting the matrix. We considered aggregates with a relative 
ice mass fraction Mi ce / (Mice+AIdirt) from to 48 % (where 
dirt refers to silicates or organics) and a porosity from to 
0.8. 
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Table 4. Dust opacities at 3.3 mm. 



Compound Size index Size range Porosity K3.3mm (m 2 kg *) 

(mm) 



Ice fraction — % Ice fraction = 22 % Ice fraction = 48 % 



Astro, silicates 


-3.0 


10" 


-4 


1 


(0. 


0.5, 


0.8) 


(8.14, 


5.42, 


4.09) 


10" 


-2 


(3.72, 


2.26, 


2.22) 


10" 


-2 


(1.95, 


1.49, 


1.53) 


10" 


-2 


Astro, silicates 


-o.u 


1 (V 
1U 


-4_ 


i n 


(U, 


u.o, 


U.o) 


( Q QQ 


O.Uo, 


0. 10) 


1 (V 
1U 


-2 




O.OO, 


o.ou ) 


1U 


-2 




z.oy, 


z.uy ) 


1U 


-2 


Astro, silicates 


-3.0 


10" 


-4_ 


-100 


(0, 


0.5, 


0.8) 


(0.88, 


1.51, 


2.81) 


10" 


-2 


(0.92, 


1.43, 


2.15) 


10" 


-2 


(0.94, 


1.32, 


1.67) 


10" 


-2 


Astro, silicates 


-3.5 


10" 


-4_ 


-1 


(0, 


0.5, 


0.8) 


(5.39, 


3.65, 


3.43) 


10" 


-2 


(2.44, 


1.75, 


2.03) 


10" 


-2 


(1.38, 


1.25, 


1.44) 


10" 


-2 


Astro, silicates 


-3.5 


10" 


-4_ 


-10 


(0, 


0.5, 


0.8) 


(4.62, 


4.98, 


5.62) 


10" 


-2 


(3.12, 


3.08, 


2.98) 


10" 


-2 


(2.23, 


2.15, 


1.91) 


10" 


-2 


Astro, silicates 


-3.5 


10" 


-4_ 


-100 


(0, 


0.5, 


0.8) 


(1.89, 


2.48, 


3.68) 


10" 


-2 


(1.49, 


1.91, 


2.43) 


10" 


-2 


(1.31, 


1.58, 


1.77) 


10" 


-2 


Astro, silicates 


-3.5 


10" 


-4_ 


-1000 


(0, 


0.5, 


0.8) 


(0.64, 


0.92, 


1.46) 


10" 


-2 


(0.55, 


0.75, 


1.14) 


10" 


-2 


(0.53, 


0.70, 


0.98) 


10" 


-2 


Astro, silicates 


-4.0 


10" 


-4_ 


-1 


(0, 


0.5, 


0.8) 


(1.78, 


1.71, 


2.68) 


10" 


-2 


(0.98, 


1.18, 


1.79) 


10" 


-2 


(0.73, 


0.98, 


1.33) 


10" 


-2 


Astro, silicates 


-4.0 


10" 


-4_ 


-10 


(0, 


0.5, 


0.8) 


(2.48, 


2.63, 


3.50) 


10" 


-2 


(1.58, 


1.72, 


2.11) 


10" 


-2 


(1.16, 


1.30, 


1.48) 


10" 


-2 


Astro, silicates 


-4.0 


10" 


-4_ 


-100 


(0. 


0.5, 


0.8) 


(2.17, 


2.46, 


3.45) 


10" 


-2 


(1.40, 


1.69, 


2.15) 


10" 


-2 


(1.14, 


1.33, 


1.53) 


10" 


-2 


Silicate mixt. 


-3.5 


10" 


-4_ 


-1 


(0, 


0.5, 


0.8) 


(0.80, 


0.52, 


0.60) 


10" 


-2 


(0.84, 


0.68, 


0.79) 


10" 


-2 


(0.87, 


0.85, 


0.96) 


10" 


-2 


Silicate mixt. 


-3.5 


10" 


-4_ 


-10 


(0. 


0.5, 


0.8) 


(1.13, 


1.00, 


0.93) 


10" 


-2 


(1.07, 


1.02, 


0.88) 


10" 


-2 


(1.16, 


1.04, 


0.91) 


10" 


-2 


Silicate mixt. 


-3.5 


10" 


-4_ 


-100 


(0, 


0.5, 


0.8) 


(0.61, 


0.72, 


0.83) 


10" 


-2 


(0.68, 


0.77, 


0.85) 


10" 


-2 


(0.82, 


0.88, 


0.89) 


10" 


-2 


Organics 


-3.5 


10" 


-4_ 


-100 


(0, 


0.5, 


0.8) 


(0.89, 


0.97, 


1.05) 


10" 


-2 


(0.87, 


0.96, 


0.99) 


10" 


-2 


(0.90, 


0.97, 


0.99) 


10" 


-2 



If the dust properties (size distribution and chemical 
properties) do not vary within the instrumental field of 
view, the flux density (W m~ 2 Hz -1 or Jy) may be written 
as (e.g., Jewitt & Luu 1990) : 



5' a 



2k 
A 2 A 2 



r(a)Q ab s7ra 2 n(a)da, 



(1) 



where n(a) oc a q is the size distribution, q is the size in- 
dex in the coma, and a m i n and a max are the minimum 
and maximum grain radii. T(a) is the temperature of the 
grains which, for the purpose of modelling thermal emis- 
sion at 3.3 mm, can be assumed in first ap proximation to 
be independent of size (| Jewitt fc Luulll990l ) and described 
by the blackbody equilibrium temperature of fast-rotating 
bodies of 174 K at 2.45 AU. The Ray leigh- Jeans approxi- 
mation for the blackbody grain emission applies. The flux 
density is then related to the mass of emitting grains M 
through the so-called dust opacity k\: 



S\ — 
with 

n\ = 



2kTMn x 
A 2 A 2 : 



XwT Qnb s na 2 n(a)da 



(2) 



(3) 



where p is the density of the grains. Densities of 1000, 1500, 
2500 kg m -3 were taken for ice, organics, and silicates, re- 
spectively. To compute k\ , we used the effective density of 
the grains, which depends on the relative proportions of the 
materials and their porosity. For example, p = 723 kg m~ 3 
for silicate grains with 50% ice content and a porosity of 
0.5 (nominal composition considered in Sects. 4.2 and 4.3). 
Sample opacities at A = 1 mm for non-porous grains made 
of various (one compon ent) materials have been presented 
bv I Jewitt fc Luul (|l992). Table H presents dust opacities at 
3.3 mm (fc3.3 mm ) for different ice/silicate mixtures, porosi- 
ties and size distributions. The maximum size considered 
in the Table is a max = 1000 mm, as evidence was found 



for cm -sized particles and even large-size fragm ents in 17P 
coma (|Reach et alJl2010t IStevenson et al.ll201Q[ ). When ice 
is included in the aggregates, the opacity ranges from 0.005 
and 0.037 m 2 kg -1 in the range of considered size indexes. 
For a max > A = 3.3 mm, K3.3 mm decreases with increas- 
ing a max because the larger particles, which are the most 
efficient radiators at millimetre wavelengths, contribute to 
the mass faster than they contribute to the radiating cross- 
section (Fig. [5]). On the other hand, K 3 . 3mm shows little 
dependence with a m i n . The absorption efficiency shows res- 
onances at particle radii between 1 and 10 mm (i.e., at ~ A, 
Fig.©, which are smeared out when the porosity increases. 
When the size distribution encompasses these resonances 
{imax = 100 and 1000 mm in the Table), porosity enhances 
the dust opacity at 3.3 mm when the "hollow sphere" ap- 
proximation is used (Fig. [5]), whereas it does not influence 
the dust opacity when vacuum is taken to be the matrix 
in the Maxwell Garnett theory. For grains made of astro- 
nomical silicates, K3.3 mm decreases for increasing amount of 
ice, as ice is more transparent than these silicates. K 3 . 3mm 
is also smaller for grains containing organics or crystalline 
silicates (i.e., forsterite). 

Applying Eq. [5] to the flux density of ~ 2.4 mJy/beam 
observed on 27 October, we derive a dust mass M = 0.8- 
6 x 10 11 kg for dust opacities in the range 0.005-0.04 m 2 
kg -1 . Obviously, this mass is a lower limit to the dust mass 
produced by the outburst, since the beam encompasses only 
part of the coma. 

Using the same approach, we computed opacities at 
0.45, 0.8 and 1.2 mm which might be useful for the interpre- 
tation of comet observations carried out with the Atacama 
Large Millimeter and submillimeter Array (ALMA) (online 
Table 0. 

4.2. Time-dependent modelling 

The in-depth analysis of the data requires a time-dependent 
model of the dust coma. Indeed, since the outburst pro- 
duced dust particles and chunks with a broad range of ve- 
locities, the dust size distribution evolved in the coma, both 
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Table 6. One-component model: best fit for astronomical silicates with 0.5 porosity and 48% ice fraction. 



Onset 
(Oct. UT) 


Size 
index 


Q>max 

(mm) 


V 
(m s- 1 ) 


Y 2(a) 


V280ct/V270ct 

r uv = 20 m 


-^280ct 
^*270ct 


ji fbcam(fa,c) 
J "dust 

(10 11 kg) 


artot (c,d) 
JU dust 

(10 11 kg) 


✓rtot (e) 
Wust 

(10 12 m 2 ) 


zo.o 


■> 

o 


i on 


1 o 




n qi 


u.uo 


^ 7 


o.o 


n 070 — O 08*3 
U.UfU u.uoo 




-3.5 


100 


60 


1.9 


0.93 


0.73 


3.0 


5.4 


1.8 - 5.8 




-4 


100 


40 


2.4 


0.84 


0.62 


2.1 


9.1 


81 - 802 


23.8 


-3 


100 


85 


3.4 


0.90 


0.61 


3.7 


5.7 


0.068 - 0.082 




-3.5 


100 


65 


2.3 


0.88 


0.64 


3.1 


5.2 


1.7 - 5.5 




-4 


100 


35 


1.7 


0.87 


0.65 


2.3 


7.9 


71 - 708 


23.8 


-3 


1000 


200 


2.3 


0.89 


0.63 


20 


27 


0.039 - 0.045 




-3.5 


1000 


125 


1.6 


0.86 


0.64 


11 


15 


1.6 - 5.2 




-4 


1000 


40 


1.5 


0.85 


0.66 


3.5 


9.4 


71 - 696 


Data 










0.84 ± 0.11 


0.80 ± 0.05 






55^ 



a) Reduced \ 2 from the fit of 27 and 28 October data. 
b ' Dust mass within the synthesized beam on 27 October. 

c ' Models with a m i„ — 1 fim. Values obtained with a m i n = 0.1 (im differ by less than 20%. 
d ' Total dust mass produced by the outburst. 

• e ' Total scattering cross-section. Calculations are for a m i„ in the range 0.1-1 fim. 
■f' From optical observations (|Li et al.ll201lT ). 




Fig. 5. Dust opacities at 3.3 mm for grains made of as- 
tronomical silicates. Values for 0% ice and porosity = 0, 
0% ice and porosity = 50%, 48% ice and porosity = 50% 
are plotted with black, red, and blue symbols, respectively. 
Plain, open circles, and crosses are for size indexes q — -3, 
-3.5, and -4, respectively. 



Fig. 6. Product Q a bs(a) X ira 2 x n(a)da for a size distri- 
bution in a~ 3 ' 7 and absorption efficiencies at 3.3 mm of 
astronomical silicates with 0.5 porosity and 48% ice frac- 
tion. 



with time and distance to nucleus. The velocity partitioning 
resulted in enhancing the relative fraction of large particles 
with respect to smaller ones in the inner coma, and might 
be observable on the visibility curves which sample a broad 
range of distances to nucleus. Alternatively, the visibilities 
at the different baseline lengths (i.e., wv-radius) may sample 
particles produced with similar kinematic properties but at 
different times. 

The number of particles injected in the coma by the 
outburst as a function of time is modelled by a function 
G(t) peaking at t = t (23.3 October UT, according to 
Hsic h et al.ll2010[ ). G(t) is described by the combination of 
two half Gaussians peaking at equal values, with widths 
AT ou tburst for t > t , and AT out burst/4 for t < t . The ef- 
fective width of G(t) is thus 5/8 x AT out b U rst- This choice 



of function was initially motivated to accelerate the conver- 
gency of the computations. 

The size distribution of the dust injected by the out- 
burst follows a power law n(a) oc a 9 , defined by the size 
index q. The velocity of the dust particles is supposed to 
vary according to : 

VdustO) = (a/ao)- 5 , (4) 

where Vq is the velocity of a grain of radius ao — 1 mm. We 
imposed Vd us t(a) < 1 km s _1 . It is worth mentioning that, 
for steady state dust production, the dust size index in the 
coma (that used in Sect. l4~T1 for computing dust opacities) 
is equal to q + 6. 
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We assumed a spherical dust coma. The density distri- 
bution of dust grains of size a is given by : 

n dus t(a,r) = — T~r^G{t - t - — — ), (5) 

where r is the distance to the nucleus and t the time of 
the observation. The dust thermal emission by unit of solid 
angle [Jy sr _1 ] at impact parameter D from the nucleus 
is computed by integration along the line of sight (LOS) 
according to: 

F X (D) = —- / / 7ra 2 Q abs (a)ri dus t(a, r)n(a)dadz. (6) 

A J LOS Ja mirl 

Visibilities V as a function of itv-radius r uv are then 
computed from this brightness distribution, taking into ac- 
count the fi nite size (with g aussian shape) of the primary 
beam (see iBoissier et al.l l2007t iBockelee-Morvan et al.l 
l2009h . The optical properties of the individual grains 
are supposed not to vary with time, though sublimation- 
induced variations might be expected. The absorption effi- 
ciencies Qabs(o) are computed as explained in Sect. 14.11 

In our calculations, the number of bins in particle sizes 
is typically 200-500. A high number of bins is a requisite 
for getting a continuous dust spatial distribution when sum- 
ming Eq. [5] over size. Actually, the number of bins is fixed 
by the condition: 

(V dust (i + 1) - V dust {i))/V dnst (i) « AT outburst /(t - *o)(7) 

to be fulfilled for all bins i. 

The velocity of the dust particles is a critical param- 
eter. Information is available for the small particles, from 
the evolution of the structure morphology of the expand- 
ing shell observed in optical images. The expansion of the 
edge of the shell is consistent with a velocity of ~ 550 
m s -1 pertaining to particles in the 0.1-100 fim range 
(|Lin et al.ll200llHsieh et al.ll2010D . However, little is known 
on the kinematics of large grains. Particles are likely ac- 
celerated by their interaction with the flow of gas liber- 
ated by sublimating icy grains. For steady-state comas, 
the terminal velocity sca les as a~ 0,5 for intermediate sizes 
(jCrifo fc R odiouov 1997). In order to roughly estimate the 
range of expected velocities Vq for 1-mm sized particles we 
performed a 1-D dust dynamics modelling specific to comet 
17P/Holmes using gas-drag coefficient for free molecular 
flow. Nucleus gravity is computed assuming a nucleus den- 
sity of 1000 kg m -3 . Gas dynamics is treated separately 
and assumes isentropic expansion. The gas total produc- 
tion rate was taken equal to 5 x 1 30 s -1 , based on the 
peak HCN production measured by iBiver et al.1 (|2008l) on 
25.9 October UT and the HCN/H 2 Q r elative abundance 
determined by iDello Russo et all ((2008) . Under these as- 
sumptions, the terminal velocity of 1-mm radius particles 
is found to be 70 m s _1 , with a size dependence for a = 100 
/ttm to 10 cm close to 5 — 0.5. The velocity of 1-100 /itm 
particles ranges from 450 to 720 m s _1 , approaching the gas 
velocity (800 m s _1 ) for the smallest sizes. The dynamics 
of 17P/Holmes's coma after the explosive event was obvi- 
ously more complex than this simple steady-state model. 
Because gases outflow more rapidly than grains, the local 
gas pressure in the dust ejecta cloud might have fastly de- 
clined with time, after the peak of gas production, reducing 
dust acceleration. For a gas production 10 times lower, the 



steady-state model predicts a terminal velocity for 1-mm 
particles ~ 3 times lower. For the calculations, we have as- 
sumed 5 — 0.5 and considered velocity values Vq between 
10 and 250 m s _1 . 



5. Results 

5.1. One-component models 

Simulations were performed with an outburst of short du- 
ration (AT out burst = 0.3 d) compared to the elapsed time 
between outburst onset and PdBI observations (see fur- 
ther discussion at the end of Sect l5.ll) . We investigated 
several outburst onset times: 23.3 UT, 23.8 UT and 24.3 
October UT. Ind eed, though the ons et time was established 
to be 23.3 UT bv lHsieh et al.l (|2010t ). this is somewhat dis- 
puted (Z. Sekanina argues for an outburst oc curring 0.4 d 
later, see note added in proof in lLi et al"1l201lD . In addition, 
the maximum rate of bright ening at optica l wavelenths oc- 
curred on 24.5 October UT (|Li et al.ll20li . 

The simulations were performed for astronomical sili- 
cates mixed with water ice (48% in mass) and a poros- 
ity of 0.5. Figures [7k.-b plot, as a function of the velocity 
parameter Vq and for S — 0.5 (Eq. H]), several quantities 
characterizing the evolution of visibilities from 27 to 28 
October: the ratios V270ct/V280ct at r uv — 20 m and 70 
m, and the ratio of the flux density within the interfero- 
metric beam F280ct/-F270ct- Also plotted is the slope index 
q v for 27 October, defined according to V270ct °e r~£ v . The 
observed values of these quantities are also given in the 
figure. We see that there is a correlation between q v and 
the decrease of F and V from 27 to 28 October. Larger- 
grain velocities lead to a faster decrease of the flux den- 
sity and to a more extended coma, which, in the Fourier 
plane, results in a faster decline of V with r uv (i.e., large 
Qv, Fig. [7J). Alternatively, for low velocities, the dust coma 
remains unresolved for most baselines and q v approaches 
the zero value characterizing an unresolved source. For a 
large enough dust parameter Voj the material contributing 
to the 3.3 mm emission forms a ring of particles which inner 
boundary has a size comparable or larger than the interfer- 
ometric beam, and the real part of the visibility becomes 
negative at large r uv . This is unlike the measurements, im- 
plying that low- velocity particles present near the nucleus 
are detected with the longest baselines at the two dates. 
Model results shown in Fig. [7] consider maximum size par- 
ticles a max equal to 100 and 1000 mm. When a max is set to 
10 mm, only a small size range is contributing to the 3.3 mm 
radiation (Fig. [6]). This thin shell of radiating grains rapidly 
moves through the field of view. Visibility curves then dis- 
play either a small r uv dependence (unresolved shell, low 
q v ) or the shape characteristics of a ring-like structure ex- 
plained above (steep r uv dependence and real part with 
negative values at large r uv ). For to — 23.8 h, the transi- 
tion between the two cases is obtained for Vq ~ 20-30 m 
s _1 (consistent with Vo X (i — to) ~ FOV radius). 

The observed visibilities were compared to the modelled 
curves. The resulting reduced xl f° r the 27 and 28 October 
data, and for the fit of both 27 and 28 October data, are 
shown in Fig. [7fc-d for to — 23.8 h and a max = 100 and 
1000 mm. The velocities Vq providing the best xt are given 
in Table [6] Vq is lower for a steeper size distribution be- 
cause the 3.3 mm emission is then sampling a relatively 
larger amount of small dust particules with larger veloci- 
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Fig. 7. Model results as a function of dust velocity parameter Vq for astronomical silicates with 0.5 porosity and 48% 
ice fraction. The outburst onset is t = 23.8 UT and the maximum dust size is a max — 100 mm (top figures (a) and (c)) 
and a max = 1000 mm (bottom figures (b) and (d)). Results for size indexes of q — -3 and q — -4 are plotted respectively 
with solid, and dotted lines, and full circles, and crosses. Left: visibility curves and interferometric-flux characteristics: 
V270ct/V280ct at r uv = 20 m (red) and 70 m (blue); interferometric flux ratio -F270ct/^280ct (green); index q v of visibility 
curve on 27 October (black). The measured values and their error bars are shown by dotted horizontal lines, and vertical 
lines, respectively. Right: reduced x 2 from data fit: 27 October data (red), 28 October data (blue), 27 and 28 October 
data (black). 




V [m s ] V [m s ] 



ties. For the same reason, Vq decreases with decreasing a max 
(Table E]). On the other hand, retrieved Vq are not very 
sensitive to the outburst onset time for to within 23.3-24.3 
UT. For example, for a max = 100 mm, q = -3/-3.5/-4, , 
V is ~ 75/60/40, 85/65/35, and 100/85/40 m s" 1 for t = 
23.3, 23.8, and 24.3 October UT respectively. Altogether, 
Vq values range from typically ~ 40 to 100 m s -1 , with 
values as large as 200 m s _1 being obtained for the size 
distribution in a~ 3 (Table [6]). Such high values are likely 
unrealistic as the velocities for 1-mm grains would then be 
comparable to the yelocities measured for ^m-sized g rains 
from optical data (jLin et alJ[200l iHsieh et all 120101 ) . As 
indicated by the reduced xti the observations are best fit- 
ted for steep size distributions (Table [B]), with Vq ~ 30-40 
m s . Figure [8] shows an example of model fits for a size 
distribution in a~ 4 . 

While satisfactory fits are obtained for the visibility 
curve obtained on 27 October (xl ~ 1-2-1.5), the fit for 28 
October is poorer (Figs. [5J . In fact, this one-component 
model fails in explaining the small temporal evolution of 
V(70 m) (see Fig. [8]) and the measured ratio F 2 soct /^270ct 
of 0.80 ± 0.05 (see Table [5]) ■ Since replenishment of the 
coma (e.g., by fragmentation processes) proceeded after 



outburst onset, we coul d question the p roduction function 
G(t) used in our model. iLi et all (|201lfl observed that the 
rate of change of the scattering cross-section followed a 
Gaussian-like function of ~ 0.4 d width and peaking on 
24.54 October UT. Our G(t) function has the same width 
(note that model results are not much sensitive to the dura- 
tion of injection of material, provided this duration remains 
small with respect to the elapsed time between outburst 
onset and the measurements). Model calculations with to 
~ 24.5 October UT do not explain the measurements. In 
other words, the production of the 3-mm radiating parti- 
cles did not proceed as for the micron-sized particles. The 
longest baselines probe typically distances from the nucleus 
r < 4000 km. The constancy of V for r uv > 60 m over Oct. 
27-28 suggests that the outburst produced a second pop- 
ulation of very slowly moving grains which 3-mm emission 
remained unresolved by the interferometric beam on both 
dates. This population of slowly moving grains produced by 
the outburst has been identified in m id-infrared images and 
referred as to the "core" component (Re ach et al.ll2010l ). Its 
origin will be discussed in Sect. |5J In the next section, our 
model will incorporate this second population of debris. 
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Table 7. Two-component model: best fit for astronomical silicates with 0.5 porosity and 48% ice fraction. 



Onset 

(October UT) 



Shell 



Size 
index 



(mm) 



V 

(m 



(10 



bcam(a,b) 
dust 
11 



kg) 



(10 



tot(b,c) 
dust 
11 



kg) 



(10 



tot(d) 
dust 
12 „2\ 



Core 
3-mm Flux 
(mJy) 



Y 2 (e) 



23.3 
23.8 



-3.5 


100 


125 


1.6 


-4.0 


100 


75 


1.4 


-3.5 


10 


60 


0.60 


-4.0 


10 


50 


0.70 


-3.5 


100 


125 


1.7 


-4.0 


100 


85 


1.4 


-3.5 


1000 


250 


8.3 


-4.0 


1000 


75 


2.9 



5.5 
9.8 
3.2 
8.4 
5.0 
9.7 
16 
9.7 



1.8 - 5.8 
89 - 873 
3.3 - 10 
95 - 924 

1.6 - 5.2 
87 - 861 

1.7 - 5.4 
70 - 715 



1.41 
1.26 
1.60 
1.57 
1.27 
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a ' Dust mass within the synthesized beam on 27 October. 



For a„ 



1 /im. Values obtained with a m i n = 0.1 /im differ by less than 20%. 



c ' Total dust mass in the shell produced by the outburst. 
d ' Total scattering cross-section for a m i„ in the range 0.1-1 /jm. 
e) Reduced y 2 fr° m the fit of 27 and 28 October data. 




uv-radius [m] 



Fig. 8. Visibility curves with the one-component model, for 
an outburst with AT out burst = 25000 s, a size distribution 
in a~ 4 with a max — 100 mm, a velocity law Vdust(fl) = 
Vq (a/ao) -0 ' 5 with Vq = 35 m s _1 being the velocity of 1- 
mm particles, and an outburst onset to = 23.8 UT. Data 
are shown by filled circles with error bars. The solid curves 
show the model results. Black (red) symbols and curves are 
for 27 and 28 October, respectively. Dust particles are made 
of a porous (50% porosity) mixture of astronomical silicates 
and ice, with an ice mass fraction of 48%. 

5.2. Two-component models 

In order to characterize the unresolved population of slowly 
moving grains, we fitted the observations by the combina- 
tion of emission of outflowing dust grains (hereafter ref- 
ered as "shell"), and a point-source emission (the "core" 
component). The shell emission is modeled as described 
in Sects. 14.21 and 15.11 The flux density from the core was 
determined from the constraint that the core remained un- 
resolved on 27 and 28 October. Figures [9^.-0 show, for two 
sets of parameters, the modelled visibilities for the two com- 
ponents and the comparison to the measurements. At r uv 
> 60 m, the ratio V270ct/V280ct is 1-1.2, consistent with 
the measured value of 1.05 ± 0.15 (see Fig. [7]). The best 
fits are obtained for grain velocities in the shell higher than 



determined from the one component model. The results of 
this two-component modelling are summarized in Table [JJ 
The flux density from the core is in the range 0.9-1.6 mJy. 
The highest values are obtained when a max = 10 mm. In 
this case, the 3.3 mm-radiating shell has a small thickness 
and its inner boundary is larger than the interferometric 
beam, so that most of the emission detected with the long 
baselines is from the material composing the core (Fig. 
Residuals between modelled and observed maps are shown 
in Fig. [TU] for the set of model parameters of Fig. [5b- 



The PdBI observations do not provide strong con- 
straints on the kinematic properties of the dust particles 
composing the core. As already discussed in Sect. [3l the 
mean velocity Vj^™ of the core should not exceed 10 m s _1 
to remain unresolved by the PdBI on the two dates. Instead, 
we used published continuum measurements obtained at 
250 GHz (i.e., 1.1 mm wavelength) using the IRAM 30-m 
telescope (HPBW = 1 1") which cover the 16 November 
to 18 December period (jAltenhoff et alJ[2009|) . Absorption 
coefficients were computed at 1.1 mm, following Sect. 14. ll 
Figure [11] shows the measurements of the dust continuum 
flux at 1.1 mm, where the 3.3 mm flux densities recorded at 
the PdBI were extrapolated to 1.1 mm, based on our cal- 
culations of dust opacities at the two wavelengths. Because 
they were performed late after the outburst, the IRAM 30- 
m telescope detected mainly emission from the low-velocity 
particles composing the core (see the modelled time evolu- 
tion of the shell emission in Fig. ITTj) . This is in agre ement 
with the conclusion obtained by iReach et al.l (|2010[ ) com- 
bining IRAM and Spitzer data sets. The time evolution 
of the core emission depends on the particle size distribu- 
tion and velocity field. Examples of model fits for q = -3.2 
(Sect. I5~3"]) are given in Fig.QTJ The velocity of 1-mm parti- 
cles in the core is typically Vq = 7-9 m s . The velocity of 
the largest particles (0.1 m) is on the order of the escape ve- 
locity (estimated to 0.8 m s _1 ). Fro m mid-infrared ima ges 
of the core obtained in March 2008. lReach et al.l ()2010t ) in- 
ferred a velocity of 9 m s _1 . The agreement with our result 
is excellent. 
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uv-radius [m] uv-radius [m] 



Fig. 9. Visibility curves for the two-component model. Model curves from the shell emission and residuals (data minus 
shell emission) are shown by a green (blue) line and green (blue) plain circles, respectively, for 27 (28) October UT. 
Parameters for the shell emission are: outburst with AT out burst = 25000 s, <o = 23.8 UT, and a size distribution in a~ 4 : 
(a) Vq = 50 m s _1 , a max = 10 mm; (b) Vq = 85 m s _1 , a max — 100 mm. The point source (core) emission is shown 
with the dashed black line. The sum of the shell and core emissions are shown in solid lines. Data are shown by filled 
circles with error bars. Black (red) symbols and curves are for 27 and 28 October, respectively. Dust particles are made 
of a porous (50% porosity) mixture of astronomical silicates and ice, with an ice mass fraction of 48%. 



5.3. Dust masses and size indexes 

The time-dependent model allows us to compute not only 
the mass of dust particles within the field of view at 
the time of the measurements M^*™, but also the total 
mass injected by the outburst M^ st . These quantities are 
given in Tables [6] and [7] for the one-component and two- 
component models, respectively. For the shell, M^°* t differs 
from M^^ n measured on 27 October by a factor of 2 to 12, 
the largest difference being obtained for the size distribu- 
tions with size index q = — 4 (Table [7]) . Most of the initial 
mass is then comprised in small, rapidly moving particles 
which are outside the PdBI field of view on 27 October. 
The dust mass found in the shell component ranges from 
3 to 16 x 10 11 kg (Table ED. However, models requiring ve- 
locities as large as 250 m s _1 for 1-mm particles can be 
excluded (Sect. 14. 2(1 . Therefore, values in the range 3-10 
x 10 11 kg are more likely, under the assumption that the 
optical properties assumed for our calculations are repre- 
sentative of comet grains. 



20 - 27 Oct. B0Q7- 



P> 



- 



-20 - 



O 



Or 



_l I I I l— 




20 -20 

RA ["] 

Fig. 10. Residuals obtained from the difference between 
observations and the two-component model. Model param- 
eters are those of Fig. [5b. Left: residuals for 27 October. 
Right: residuals for 28 October. The levels are la spacing 
with cr = 0.093 and 0.12 mJy/beam for October 27 and 
28, respectively. The synthesized interferometric beam is 
plotted in the left corner. 



For comparison with optical data, Tables H2 and [7] 
present total dust scattering cross-sections C^°* t inferred 
for the different models, assuming minimum particle sizes 
in the range <z m j n = 0.1-1 /im. Particles with sizes less 
than 0.1 /im are not considered as they do not scatter effi- 
ciently visible light. The scattering cross-section measured 
from optical observations reached 5.5 x 10 13 m 2 at the peak 
of brightening on ~ 25.0 October UT (|Li et al.ll20TTh . For 
o-min = 1 A 4111 ) the best fit for the two-component model 
implies C^ st = 1.6-3.3 x 10 12 m 2 for q = -3.5, and C£f st 
= 70-95 x 10 12 m 2 for q = -4 (Table EJ), which suggests a 
size distribution in the shell component close to a~ 4 if the 
minimum size is 1 /im. From interpolation, we found that 
the size index fitting both optical and radio data is q ~ -3.7 
for a m i n — 0.1 /im and -3.9 for a m in = 1 fni. These values 
are obtained for astronomical silicates mixed with water 
ice and a porosity of 50%. Taking into account the range of 
grain compositions and porosities considered in Tabled q 
is within [-3.8, -3.6] and [-4.1, -3.8], for a m i n — 0.1 and 1 
/im, respectively. The inferred size index is not significantly 
sensitive to a max . However, our model is not considering 
time-dependent dust production processes, such as particle 
fragmentation, responsible for the rapid i ncrease of the scat- 
tering cross-se ction after outburst onset (Hsich e t alj|20iot 
iLi et alJl20ll . This introd uces additional uncertainties in 
the size index determination. 

The size distribution in the core component can be eval- 
uated from optical data contemporaneous to PdBI data. 
We assumed that the bright condensation surrounding the 
nucleus seen in optical images is the visible counterpart of 
this core component. From images obtained by N. Bivcr 
and J. Nicolas (private communication) on 26.9 and 27.8 
October UT, the scattering cross-section in a 6" diame- 
ter aperture corresponding to the PdBI field of view was 
0.7 and 0.4% the total scattering cross-section, implying 
a scattering cross-section of at most ~ 3.5 x 10 11 m 2 in 
the core on 27.1 October UT. The small variation of the 
cross-section from 27 to 28 October can be interpreted as 
indicative of slowly moving particles (~ 15 m s^ 1 , if we as- 
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sume that dust production from the nucleus or from frag- 
mentation was negligible between the two dates). The Q a bs- 
weighted cross-section of 3 mm-radiating particles derived 
from the 1.3 mJy flux density (Table© is 1.8 x 10 9 m 2 . The 
comparison with the scattering cross-section shows that the 
core contained a larger proportion of large particles, com- 
pared to the shell. Assuming a m i n — 0.1 /mi and a max 
1 to 1000 mm, optical and radio data can be reconciled 
with a size index of -3.2 to -3.0. The results are not much 
sensitive to the minimum particle size: q is within [-3.3, 
-3.0] for a m i n — 1 /im. Using these size indexes, we derive 
dust masses for the core of 0.7, 1.2, and 3.7 (x 10 11 ) kg for 
a-max — 10, 100, and 1000 mm, respectively. Optical and ra- 
dio data cannot be both explained for maximum grain sizes 
set to < 2 mm. These calculations pertain to the nominal 
composition (50% astronomical silicate, 48% ice, porosity of 
0.5). Considering the range of dust composition and poros- 
ity in Tabled q remains within [-3.35, -3.0] for a m i n = 0.1 
/im. Given uncertainties in the size distributions for the two 
components, especially on a max , the ratio A^ust /Mdust' i s 
in the range 0.1-1. 

In summary, we derived dust masses of 3-10 x 10 11 
kg and 0.7-4 x 10 11 kg for the shell and core compo- 
nents, respectively. These calculations are valid for astro- 
nomical silicates mixed with water ice, with a porosity of 
0.5. Additional uncertainties come from the optical proper- 
ties of the dust ejecta which are not precisely known. Based 
on our calculations of dust opacities for different materials 
(Table S]), dust masses could be 60% higher (material re- 
ferred as to organics and silicate mixture in the table) or 
up to a factor 2.6 lower (higher porosity, ice- free). 

6. Discussion 

We performed a detailed analysis of the 3.3 mm dust ther- 
mal emission from comet 17P/Holmes observed 4-5 days af- 
ter its massive outburst with the IRAM Plateau de Bure in- 
terferometer. The data were combined to optical and IRAM 
30-m telescope measurements to characterize the properties 
of the ejecta cloud, namely the dust size distribution, its 
kinematics and mass. 

Two distinct dust components, with different kinematic 
properties and size distributions, arc identified in the data. 
The large- velocity component, with typical velocities Vo of 
50-100 m s^ 1 for 1 mm particles, corresponds to the fast 
(550 m s _1 ) expanding shell observed in optical images and 
displays a steep size distribution with a size index estimated 
to q = -3.7. The size dependence of the velocity is consis- 
tent with gas drag fSect. l4T2"]) . The very high gas production 
rates measured shortly after the outburst (> 10 30 s _1 ) and 
their fast temporal d rop-off argue for gas production from 
subliming icy gr ains (jBiver et al.l 120081 : ICombi et al.ll2007t 
ISchleichenl200il . High gas pressures caused the particles 
constituting the shell to be accelerated to large terminal 
velocities. The second component, the "core", consists in 
slowly-moving particles with kinematic proper ties (Vn = 7- 
9_m_s _1 ) consistent with results obtained by iReach et al.l 
l|2010t l from the study of mid-infrared images. Constraints 
obtained on the size index and maximum dust size show 
that the core is dominated by large (a max > 2 mm) parti- 
cles with a shallow s ize index \q\ < 3.4. The same conclu- 
sion was obtained bv lReach et a l. (2010) from the analysis 
of the Spitzer images. In addition, the core revealed fea- 
tureless 10- /im spectra in mid-November 2007, in contrast 




20 40 
Time from outburst onset [day] 



Fig. 11. Evolution of the 250 GHz flux f rom IRAM 30-m 
telescope (red dots. lAltenhoff et al.l 12009^ and PdBI data 
(black dots). PdBI 90-GHz fluxes were converted to 250 
GHz fluxes, applying a conversion factor of 22.1 determined 
by thermal modelling. The beam size is HPBW = 11" for 
the 30-m data and ~ 6" for the PdBI data. The observed 
core (Table [7]) and shell emissions on 27-28 October are 
shown by opened circles and triangles, respectively. The 
dotted-dashed line shows the model for the shell emission 
(q = -4, a max = 100 mm, V Q = 85 m s" 1 ) with HPBW = 
6" for the two first dates, and 11" later on. Two-component 
models (shell+core) are shown with solid, dashed and dot- 
ted lines, with the core parameters being q = -3.2 and 
(a max = 100 mm, Vo = 5 m s -1 ), (a rnax = 100 mm, V a = 9 
m s _1 ), and (a max — 10 mm, Vq — 7 m s _1 ), respectively. 



to the sh ell, consistent with much larger grains populating 
the core (|Reach et al.ll2010ft . 

The origin of the low-velocity core component may be 
related to the fast decline of the gas pressure with time 
in the inner coma resulting from the fast (> 1 km s _1 , 
iBoissier et al.l I2009T) expansion of the first generation of 
gases. Under such conditions, we expect gas velocities to 
vary both with time and distance to nucleus, with the older 
molecules present at larger distances re aching higher ve- 
locities, and this w as indeed observed (jBiver et al.l 120081 : 
IBoissier et al.l [20091) . The fast decline of gas pressure and 
velocity in the inner coma resulted in a strong decrease of 
the gas aerodynamic force acting on grains, making their 
acceleration less effective. Grains produced after the peak 
of gas production were thus accelerated to lower terminal 
velocities. Another effect to consider is that efficient grain 
acceleration requires gas pressure gradient. Particles in the 
centre of the expanding cloud were subject to limited accel- 
eration since the gas pressure gradient caused by the sub- 
liming grains was minimum in this region. Finally, slowly- 
moving particles, in the trailing side of the ejecta cloud, 
were decelerated due to the positive gradient of pressure 
in the radial direction combined with nucleus gravity. The 
kinematics of the most massive debris was possibly not sig- 
nificantly affected by gas-drag; providing they were not ac- 
celerated by rocket forces, their velocity remained close to 
their separation velocity from the nucleus. Thus, the low ve- 
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locity of the particles and debris composing the core might 
have different origins. 

We derived dust masses of 3-10 x 10 11 kg and 0.7-4 x 
10 11 kg for the shell and core, respectively, using grains 
made of astronomical silicates mixed with water ice as 
cometary analogues. The ratio M^/M^ is in the range 
0.1-1, lower than the value ~ 2 determined by iReach et al.l 
(2010J). The large mass in the shell component and its steep 
size distribution indicate that the disruption of the nucleus 
was accompanied or followed by a massive production of 
small particles, indicative of material with small cohesive 
strength. 

We found that the total dust mass Md us t injected by 
the outburst was in the range 4-14 x 10 11 kg. This is in 
the high end of reported values based on optical or infra red 
observations. From optical observations. iLi et al. ( 2011) es- 
tim ated M dust to 0.2 -9 x 10 11 kg , while ISchleicherl (|2009) 
and ISekaninal (120081) favour val u es of 1-2 x 10 11 kg. The 
estimate made by IReach et al.l (|2010l) from their Spitzer 
thermal observations is 0.1 x 10 11 kg. The inconsistency 
between these concurrent measurements reflects different 
assumptions on the size distribution, especially on the max- 
imum particle size in the ejecta. For example, IReach et al.l 
(2010!) made their estimations using characteristic sizes in 
the core, blob and shell components of 200, 8 and 2 /zm, re- 
spectively, so that their values are lower limits to the dust 
masses. 

The total dust mass is about 10 times higher than the 
integrat ed water production of ~ 5 x 10* ° kg over several 
months (jBiver et al.ll2008tlSchleicherll2009D . Most of the wa- 
ter was produced on a timescale of the order of the day af- 
ter outburst onset from small dirty icy grains. O n-going gas 
produ ction was observed through March 2008 (Schleicher 
2009) . The spatial distribution of OH radicals suggests out- 
gassing from the inne r part of the com a, typically < 1-2 x 
10 4 km according to ISchleicherl (|2009l ). A significant con- 
tribution from the icy debris in the core is likely, though 
detailed modelling is required to investigate whether this 
interpretation is consistent wi th the observed r adial distri- 
bution of OH. As discussed bv ISchleicherl (|2009f ). the source 
of the low production observed in March 2008 is more likely 
the nucleus. 

Assuming a nucleus bulk density equal to the dust 
density (700 kg m~ 3 , Sect I4.2[) . the mass ejected by the 
outburst is 3-9% the nucleus mass, and corresponds to a 
cube having a side length of 0.8-1.2 km. In other words, 
the dramatic outburst experienced by 17P/Holmes was 
caused by a massive disruption of part of its nucleus. 
Possible me c hanism s for comet splitting are discussed by 
iBoehnhardtl (|2005ft . The simple fact that 17P/Holmes also 
experienced massive outbursts in 1892 and 1893 indicates 
that the outburst was not induced by a collision with an- 
other body. The idea that the outburst was triggered by the 
exothermic pha se transiti o n of b u ried amorphous wat er ice 
is defended by ISekaninal (|2008D . IReach et all (j2010D . and 
ILi et al.l (1201 lh. However, from thermal m odelling of nu- 
cleus interior, Kossack i fc Szutowicd (|2010D reach the con- 
clusion that this process is unable to explain such a massive 
fragmentation. The origin of the fragmentation of comet 
17P/Holmes remains elusive. 



7. Summary 

The 3.3-mm continuum emission from comet 17P/Holmes 
was observed with the Plateau de Bure interferometer 4-5 
days after its dramatic 2007 outburst. The main results of 
a detailed analysis of these observations are: 

— The peak position of the brightness distribution coin- 
cides within 1.5cr with the nucleus position. Some ex- 
cess emission is detected southward the nucleus position 
on 27 October. This excess emission disappeared on 28 
October. 

— Two distinct dust components with different kinematics 
properties are identified in the data. 

— The large- velocity component, with typical velocities Vq 
of 50-100 m s -1 for 1 mm particles, displays a steep size 
distribution with a size index estimated to q = -3.7 for 
the nominal silicate-ice mixture assuming a minimum 
grain size of 0.1 ^m (q within [-3.8, -3.6] considering a 
wide range of mixtures and porosities). It corresponds 
to the fast expanding shell observed in optical images. 
Velocities are consistent with acceleration by gas pres- 
sure resulting from the sublimation of icy grains. 

— The slowly- moving component (Vq = 7-9 m s^ 1 ), re- 
ferred to as the core, has a shallower size index \q\ < 
3.4, compared to the shell. These particles were not ef- 
ficiently accelerated by gas drag. The dust mass in the 
core is in the range 0.1-1 that of the shell. 

— Using optical constants pertaining to porous grains 
(porosity of 0.5) made of astronomical silicates mixed 
with water ice (48% in mass) , the total dust mass Md us t 
injected by the outburst is estimated to 4-14 x 10 11 kg, 
corresponding to 3-9% the nucleus mass. The dramatic 
outburst experienced by 17P/Holmcs was caused by a 
massive disruption of its nucleus. 
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Table 5. Dust opacities at 0.45, 0.8, and 1.2 mm. 



Compound Size index Size range Porosity n\ (m 2 kg 1 ) 

(mm) 



Ice fraction = % Ice fraction = 22 % Ice fraction = 48 % 



A = 1.2 mm 



Astro, silicates 


-3.0 


10" 


-4 


-1 


(0, 


0.5, 


0.8) 


(3.99, 4.67, 


4.51) 


10" 


1 


(3.00, 2.72, 


2.23) 


10" 


1 


(2.17, 1.70, 1.43) 


10" 


\ 


Astro, silicates 


-3.0 


10" 


-4 


-10 


(0, 


0.5, 


0.8) 


(0.93, 1.58, 


2.70) 


10" 




(0.96, 1.40, 


1.96) 


10" 


\ 


(0.96, 1.23, 1.43) 


10" 


\ 


Astro, silicates 


-3.0 


10" 


-4 


100 


(0, 


0.5, 


0.8) 


(1.44, 2.73, 


5.86) 


10" 


2 


(1.70, 3.04, 


6.05) 


10" 


2 


(1.96, 3.32, 6.03) 


10" 


2 


Astro, silicates 


-3.5 


10" 


-4 


-1 


(0, 


0.5, 


0.8) 


(3.48, 3.72, 


3.70) 


10" 


1 


(2.39, 2.13, 


1.95) 


10" 


1 


(1.68, 1.39, 1.29) 


10" 


1 


Astro, silicates 


-3.5 


10" 


-4_ 


-10 


(0, 


0.5, 


0.8) 


(1.59, 2.19, 


3.11) 


10" 


1 


(1.33, 1.64, 


2.03) 


10" 


1 


(1.18, 1.33, 1.42) 


10" 


1 


Astro, silicates 


-3.5 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(0.56, 0.80, 


1.28) 


10" 


1 


(0.49, 0.67, 


1.01) 


10" 


1 


(0.47, 0.61, 0.85) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(1.56, 1.80, 


2.45) 


10" 


1 


(1.07, 1.19, 


1.53) 


10" 


1 


(0.81, 0.90, 1.09) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(1.42, 1.81, 


2.61) 


10" 


1 


(1.05, 1.28, 


1.66) 


10" 


1 


(0.87, 1.01, 1.18) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(1.22, 1.53, 


2.26) 


10" 


1 


(0.90, 1.11, 


1.49) 


10" 


1 


(0.75, 0.90, 1.10) 


10" 


1 


Silicate mixt. 


-3.0 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(1.07, 0.92, 


0.69) 


10" 


1 


(1.05, 0.81, 


0.66) 


10" 


1 


(1.10, 0.81, 0.67) 


10" 


1 


Silicate mixt. 


-3.0 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(0.42, 0.56, 


0.69) 


10" 


1 


(0.51, 0.61, 


0.70) 


10" 


1 


(0.62, 0.70, 0.73) 


10" 


1 


Silicate mixt. 


-3.0 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(0.98, 1.66, 


3.05) 


10" 


2 


(1.25, 2.08, 


3.59) 


10" 


2 


(1.58, 2.56, 4.12) 


10" 


2 


Silicate mixt. 


-3.5 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(0.84, 0.70, 


0.59) 


10" 


1 


(0.80, 0.64, 


0.59) 


10" 


1 


(0.84, 0.67, 0.62) 


10" 


1 


Silicate mixt. 


-3.5 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(0.53, 0.61, 


0.68) 


10" 


1 


(0.60, 0.64, 


0.68) 


10" 


1 


(0.70, 0.72, 0.70) 


10" 


1 


Silicate mixt. 


-3.5 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(0.24, 0.30, 


0.42) 


10" 


1 


(0.26, 0.34, 


0.46) 


10" 


1 


(0.30, 0.39, 0.51) 


10" 


1 


Silicate mixt. 


-4.0 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(0.35, 0.36, 


0.45) 


10" 


1 


(0.37, 0.40, 


0.49) 


10" 


1 


(0.42, 0.46, 0.54) 


10" 


1 


Silicate mixt. 


-4.0 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(0.37, 0.41, 


0.51) 


10" 


1 


(0.41, 0.45, 


0.54) 


10" 


1 


(0.48, 0.52, 0.58) 


10" 


1 


Silicate mixt. 


-4.0 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(0.36, 0.38, 


0.48) 


10" 


1 


(0.36, 0.42, 


0.51) 


10" 


1 


(0.42, 0.48, 0.56) 


10" 


1 


A = 0.8 mm 






































Astro, silicates 


-3.0 


10" 


_4 


-1 


(0, 


0.5, 


0.8) 


(6.04, 8.17, 


9.90) 


10" 




(5.09, 5.70, 


5.17) 


10" 


\ 


(4.19, 3.94, 3.33) 


10" 


\ 


Astro, silicates 


-3.0 


10" 


_4 


-10 


(0, 


0.5, 


0.8) 


(1.16, 2.06, 


3.93) 


10" 




(1.28, 2.05, 


3.31) 


10" 


\ 


(1.37, 1.99, 2.70) 


10" 


\ 


Astro, silicates 


-3.0 


10" 


-4 


-100 


(0, 


0.5, 


0.8) 


(1.64, 3.20, 


7.08) 


10" 


2 


(2.02, 3.68, 


7.68) 


10" 


2 


(2.38, 4.16, 8.13) 


10" 


2 


Astro, silicates 


-3.5 


10" 


_4 


-1 


(0, 


0.5, 


0.8) 


(6.02, 7.25, 


8.46) 


10" 


\ 


(4.55, 4.77, 


4.57) 


10" 


\ 


(3.56, 3.30, 3.03) 


10" 


\ 


Astro, silicates 


-3.5 


10" 


-4_ 


-10 


(0. 


0.5, 


0.8) 


(2.41, 3.47, 


5.32) 


10" 


1 


(2.14, 2.83, 


3.81) 


10" 


1 


(2.00, 2.42, 2.87) 


10" 


1 


Astro, silicates 


-3.5 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(0.81, 1.20, 


1.97) 


10" 


1 


(0.76, 1.04, 


1.61) 


10" 


1 


(0.72, 0.96, 1.41) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(3.19, 3.97, 


5.68) 


10" 


1 


(2.34, 2.79, 


3.52) 


10" 


1 


(1.89, 2.13, 2.51) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(2.74, 3.62, 


5.46) 


10" 


1 


(2.12, 2.69, 


3.59) 


10" 


1 


(1.82, 2.17, 2.61) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(2.29, 3.03, 


4.63) 


10" 


1 


(1.81, 2.28, 


3.11) 


10" 


1 


(1.54, 1.86, 2.32) 


10" 


1 


Silicate Mixt. 


-3.0 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(1.97, 1.99, 


1.64) 


10" 


1 


(2.06, 1.98, 


1.56) 


10" 


1 


(2.42, 1.95, 1.61) 


10" 


1 


Silicate Mixt. 


-3.0 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(0.63, 0.93, 


1.29) 


10" 


1 


(0.80, 1.07, 


1.38) 


10" 


1 


(0.98, 1.27, 1.50) 


10" 


1 


Silicate Mixt. 


-3.0 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(1.17, 2.12, 


4.22) 


10" 


2 


(1.56, 2.73, 


5.17) 


10" 


2 


(1.99, 3.40, 6.18) 


10" 


2 


Silicate Mixt. 


-3.5 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(1.72, 1.62, 


1.42) 


10" 


1 


(1.74, 1.60, 


1.40) 


10" 


1 


(1.98, 1.63, 1.47) 


10" 


1 


Silicate Mixt. 


-3.5 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(0.90, 1.12, 


1.36) 


10" 


1 


(1.07, 1.23, 


1.42) 


10" 


1 


(1.26, 1.41, 1.52) 


10" 


1 


Silicate Mixt. 


-3.5 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(0.35, 0.47, 


0.70) 


10" 


1 


(0.41, 0.55, 


0.78) 


10" 


1 


(0.49, 0.64, 0.89) 


10" 


1 


Silicate Mixt. 


-4.0 


10" 


-4 


-1 


(o, 


0.5, 


0.8) 


(0.81, 0.86, 


1.03) 


10" 


1 


(0.87, 0.96, 


1.13) 


10" 


1 


(1.03, 1.09, 1.24) 


10" 


1 


Silicate Mixt. 


-4.0 


10" 


-4 


-10 


(o, 


0.5, 


0.8) 


(0.76, 0.89, 


1.11) 


10" 


1 


(0.88, 1.00, 


1.20) 


10" 


1 


(1.03, 1.16, 1.31) 


10" 


1 


Silicate Mixt. 


-4.0 


10" 


-4 


-100 


(o, 


0.5, 


0.8) 


(0.68, 0.78, 


1.01) 


10" 


1 


(0.75, 0.89, 


1.09) 


10" 


1 


(0.88, 1.03, 1.21) 


10" 


1 


A = 0.45 mm 






































Astro, silicates 


-3.0 


10" 


-4 


-1 


(0, 


0.5, 


0.8) 


(0.93, 1.50, 


2.33) 






(0.94, 1.27, 


1.53) 






(0.90, 1.06, 1.05) 






Astro, silicates 


-3.0 


10" 


-4 


-10 


(0, 


0.5, 


0.8) 


(1.46, 2.74, 


5.71) 


10" 




(1.74, 3.00, 


5.67) 


10" 


\ 


(1.98, 3.23, 5.45) 


10" 


\ 


Astro, silicates 


-3.0 


10" 


-4 


-100 


(0, 


0.5, 


0.8) 


(1.95, 3.86, 


8.85) 


10" 


2 


(2.45, 4.62, 


9.99) 


10" 


2 


(2.96, 5.40, 11.2) 


10" 


2 


Astro, silicates 


-3.5 


10" 


_4 


-1 


(0, 


0.5, 


0.8) 


(1.15, 1.61, 


2.27) 






(1.02, 1.24, 


1.43) 






(0.90, 0.99, 0.99) 






Astro, silicates 


-3.5 


10" 


-4_ 


-10 


(0, 


0.5, 


0.8) 


(0.41, 0.63, 


1.03" 






1 c\ ac\ n err 

(0.40, 0.55, 


0.83" 






(0.39, 0.51, 0.69) 






Astro, silicates 


-3.5 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(1.35, 2.08, 


3.56) 


10" 


1 


(1.31, 1.89, 


3.01) 


10" 


1 


(1.30, 1.81, 2.72) 


10" 


1 


Astro, silicates 


-4.0 


10" 


-4_ 


-1 


(o, 


0.5, 


0.8) 


(0.81, 1.12, 


1.74) 






(0.67, 0.86, 


1.14) 






(0.58, 0.70, 0.82) 






Astro, silicates 


-4.0 


10" 


-4_ 


-10 


(o, 


0.5, 


0.8) 


(0.66, 0.94, 


1.51) 






(0.56, 0.75, 


1.04) 






(0.50, 0.63, 0.79) 






Astro, silicates 


-4.0 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(0.56, 0.79, 


1.27) 






(0.47, 0.63, 


0.88) 






(0.42, 0.53, 0.67) 






Silicate mixt. 


-3.0 


10" 


-4_ 


-1 


(o, 


0.5, 


0.8) 


(0.41, 0.49, 


0.52) 






(0.49, 0.54, 


0.51) 






(0.57, 0.60, 0.54) 






Silicate mixt. 


-3.0 


10" 


-4_ 


-10 


(o, 


0.5, 


0.8) 


(0.95, 1.58, 


2.69) 


10" 


1 


(1.26, 1.97, 


3.14) 


10" 


1 


(1.59, 2.42, 3.61) 


10" 


1 


Silicate mixt. 


-3.0 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(1.48, 2.77, 


5.91) 


10" 


2 


(2.02, 3.64, 


7.52) 


10" 


2 


(2.60, 4.62, 9.25) 


10" 


2 


Silicate mixt. 


-3.5 


10" 


-4_ 


-1 


(o, 


0.5, 


0.8) 


(0.41, 0.45, 


0.46) 






(0.48, 0.48, 


0.47) 






(0.54, 0.54, 0.50) 






Silicate mixt. 


-3.5 


10" 


-4_ 


-10 


(o, 


0.5, 


0.8) 


(0.18, 0.24, 


0.33) 






(0.22, 0.28, 


0.37) 






(0.27, 0.33, 0.41) 






Silicate mixt. 


-3.5 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(0.63, 0.87, 


1.35) 


10" 


1 


(0.77, 1.04, 


1.57) 


10" 


1 


(0.93, 1.25, 1.82) 


10" 


1 


Silicate mixt. 


-4.0 


10" 


-4_ 


-1 


(o, 


0.5, 


0.8) 


(2.41, 2.73, 


3.38) 


10" 


1 


(2.81, 3.18, 


3.74) 


10" 


1 


(3.31, 3.73, 4.18) 


10" 


1 


Silicate mixt. 


-4.0 


10" 


-4_ 


-10 


(o, 


0.5, 


0.8) 


(2.09, 2.53, 


3.31) 


10" 


1 


(2.50, 2.98, 


3.70) 


10" 


1 


(3.00, 3.52, 4.15) 


10" 


1 


Silicate mixt. 


-4.0 


10" 


-4_ 


-100 


(o, 


0.5, 


0.8) 


(1.79, 2.14, 


2.85) 


10" 


1 


(2.12, 2.54, 


3.21) 


10" 


1 


(2.52, 3.00, 3.62) 


10" 


1 



